Clusters under strong VUV pulses: 
A quantum-classical hybrid-description incorporating plasma effects 
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The quantum-classical hybrid-description of rare-gas clusters interacting with intense light pulses 
which we have developed is described in detail. Much emphasis is put on the treatment of screening 
electrons in the cluster which set the time scale for the evolution of the system and form the link 
between electrons strongly bound to ions and quasi-free plasma electrons in the cluster. As an 
example we discuss the dynamics of an Ari47 cluster exposed to a short VUV laser pulse of 20 eV 
photon energy. 
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I. INTRODUCTION 

The interaction of rare-gas clusters with intense laser 
pulses has attracted considerable attention over the last 
decade since spectacular experiments showed how effec- 
tive energy from the laser beam can be coupled into the 
cluster, culminating in the demonstration of fusion in 
exploding deuterium clusters Other efforts aim at 
understanding the interaction of matter beyond a few 
atoms with light of different wavelength, motivated by 
the fact that little is known about such processes with in- 
tense light of frequencies different than the infrared (IR) 
range around 800 nm wavelength. First experiments at 
FLASH in Hamburg have revealed very efficient energy 
absorption at 98 nm wavelength and have triggered 
theoretical research jS, l4Jij@l as recently reviewed Q. 

Clearly, the time-dependent description of a finite 
many-body system interacting with a laser pulse is chal- 
lenging. A full quantum description is numerically im- 
possible and would probably provide little insight con- 
ceptually. One can observe two different strategies in the 
literature to tackle this problem: (i) A more or less full 
quantum approach for a single atom in the cluster treat- 
ing the influence of the other atoms and ions in the cluster 
more approximately as a kind of environment 0, Q . (ii) 
A classical propagation of all charged particles in time 
(ions and electrons) with quantum mechanical elements 
added in form of rates on various levels of sophistication 

We have followed the latter route, describing the mo- 
tion of all charged particles classically, while they inter- 
act and are subjected to the external dipole-coupled laser 
field. Initial bound electronic motion is not treated ex- 
plicitely, but only in form of (quantum) ionization rates 
which are integrated into the classical dynamics with a 
Monte Carlo scheme, see The approach works very 
well, as long as the electrons, once ionized from their 
mother-ions, behave essentially classically. This is the 
case if they are subject to a strong external field (as it 
is the case with strong IR pulses) or are instantaneously 
free (as in the case of hard X-rays). Under VUV radia- 



tion, however, the photo-electrons stay often very close 
to their mother-ions or other ions, effectively screening 
them and modifying subsequent ionization processes. 

The reason for this behavior is the size of the quiver 
motion (typically less than a ground state electron or- 
bit of the rare-gas atom) and the small kinetic energy 
which remains of the photon energy in excess of the ion- 
ization potential. The latter is not the case for hard 
X-rays where the photon energy is high enough to re- 
move the photo-electrons completely from the cluster 
(at least for moderate cluster sizes). For intense IR- 
fields, on the other hand, the photon energy is much too 
low for the photo-electrons to leave the cluster instanta- 
neously, but the quiver motion is (for typical intensities) 
of the size of the cluster or even larger implying that the 
photo-electrons certainly will not remain in the vicinity 
of specific ions. Rather, they are dragged back and forth 
through the cluster by the laser field. Hence, the photo- 
ionization of ions surrounded by screening electrons is 
a phenomenon unique to VUV radiation, adding to the 
challenge for a theoretical description. On the other 
hand, as will become clear subsequently, exactly those 
electrons which screen individual ions define a timescale, 
suitable to formulate a coarse grained electron dynam- 
ics in the cluster. It is the key to incorporate physical 
processes which in our approach lie at the interface of 
classical and quantum description, such as the influence 
of the surrounding charged particles (ions and electrons) 
on the photo-ionization rate of an ion in the cluster. 

In the next section we summarize our quantum- 
classical hybrid description which contains the quantum 
photo-ionization rates of multi-electron ions, the classical 
propagation, and, as an element in between, the treat- 
ment and identification of those electrons which screen 
the ions. In Sect. IIIII we explain how to deal with photo- 
ionization of bound electrons into the plasma of cluster 
electrons. Section llVI discusses as an application and il- 
lustrative example the illumination of Ari47 with a VUV 
pulse of 62 nm wavelength. The paper ends with a short 
summary in Sect. [V] 
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II. HYBRID-DESCRIPTION OF 
TIME-DEPENDENT CLUSTER EVOLUTION 

The interaction of a cluster with intense radiation can 
be partitioned into three parts: (A) atomic ionization, 
(B) cooperative interaction of ions among each other 
and with electrons, and (C) relaxation. During phase 
(A) which lasts approximately until every second atom 
in the cluster is (singly) ionized, one can simply ap- 
ply the atomic photo ionization cross sections for many- 
electron atoms supplied, e.g., in ^17, 18]. Gradually, the 
bound electrons feel the cluster environment, which has 
roughly three effects (i) the ionization potential is low- 
ered through close by ions, (ii) previously ionized elec- 
trons trapped by the cluster screen the ion under consid- 
eration, (iii) the global field, generated by electrons and 
ions, modifies the ionization potential as well. We will 
treat all three effects which happen during phase (B) of 
the dynamics in Sect. Hill Here, we describe briefly how 
to extrapolate the known photo-ionization cross sections 
below the ionization threshold. 



A. Atomic photo-ionization 

The calculation of photo-ionization and excitation for 
individual energy levels of isolated atoms and ions is 
straight forward. Including the (quantum) photo rates 
into the classical propagation of charged particles in the 
cluster, we do not resolve the angular dependences of the 
photo cross section. Averaging over angular degrees of 
freedom considerably simplifies the cross sections which 
only depend on the radial wavefunctions for the respec- 
tive mean initial and final energies of the photo transi- 
tion, for a similar philosophy see 0- We start from the 
dipole matrix element for linearly polarized light along 
the z direction and consider the transition between initial 
and final atomic states with well defined orbital quantum 
numbers Im 



{f\z\i)^dr{tf,li) / dnYl^^^{n) COS 6 Y,^„.M (1) 



with the radial dipole matrix element 



drug J (r) r ue^ (r) . 



(2) 



Within the independent particle picture all states in a 
shell defined by orbital angular momentum £i are degen- 
erate and have the same radial wavefunction ue^ (r) . It 
is an eigenfunction to a spherically symmetric effective 
single particle potential obtained (along with the eigen- 
function) from the Cowan code (i7[. Within this mean 
field approximation the full photo cross section is easily 
obtained by summing over all available final states and 
averaging over the initial states in the shell £, 

agiu;) = ^^^[{i+l)dlie,i+l) + idlie,e-l)]. (3) 



Here, Wi gives the number of available initial electrons. 

Furthermore, in the cluster the ions are surrounded by 
other electrons and ions. The latter lower the potential 
barriers for ionization into the cluster. Hence, what is 
discrete photo excitation in an isolated atom becomes 
at the same photon energy ionization into the cluster. 
Therefore, we need an interpolation of the discrete exci- 
tation spectrum distributing the oscillator strength over 
the respective photon energy interval. 

In the following, we compare two different interpola- 
tions. The first one is taken from [T^ . There, approx- 
imate analytical expressions for the photo-absorption of 
many electron atoms are derived. The corresponding re- 
sults for the Ar atom are shown in Fig. [1] with dashed 
lines. For comparison, the cross sections with hydro- 
genic wavefunctions [2C| are also shown. In the second 
approximation we define a continuous photo-absorption 
in the region of discrete spectral lines by demanding 
that the renormalized photo excitation cross section 
merges smoothly with the photo-ionization cross section 
at threshold [21j. This is achieved by distributing the 
the oscillator strength /„ of a spectral line En over an 
interval half way to each of the adjacent lines such that 



cr„{uj) 



2nafn 



(En+l — En-l)/2 ' 



(4) 



where cr„(a') is now the interpolated photo-absorption 
cross section Eq. ^ for {En-i + En)/2 < uj < {En + 
En+l)/ 2. The result (solid line in Fig. [T]) shows rea- 
sonable agreement with the analytical approximation for 
low and high energies. For intermediate energies, the well 
known Cooper minima lead to considerably lower values 
than in the analytical approximation which does not ac- 
count for this interference effect. Hence, we will use in 
the following the approximation Eq. Q. 



B. Classical propagation for the Coulomb 
interaction under VUV radiation 



The propagation of the classical particles is in prin- 
ciple straight forward. More refined methods, such as 
tree-codes are only worth the effort of coding for large 
clusters (10^ electrons and ions and more). Another is- 
sue is the Coulomb interaction. Using the real Coulomb 
potential with its singularity is numerically very costly 
(small time steps close to the singularity) and leads for 
more than two bound electrons to artificial autoioniza- 
tion since one electron can fall into the nucleus (below the 
quantum ground state energy) and another one can be 
ionized with the released energy. In strong field physics 
(at IR frequencies) the so called soft-core potential 



U{r) - - 



(r2 + a2)i/2 



(5) 



has been used routinely where the singularity is cut 
off by the smoothing parameter a chosen to get poten- 
tial depths (slightly) below the true ionization potential 



3 









FIG. 1: Atomic photo-ionization rates extrapolated below the 
threshold according to [l^ (dashed) and by interpolation of 
the discrete spectrum (solid). The dashed-dotted line is an 
approximation with hydrogenic wavefunctions, see text. The 
initial configuration is displayed in each panel. 



of the atom or ion. As long as the quiver amplitude 



•^quiv 



— FjiJ^ ;2> a, the cut off is irrelevant for the dy- 
namics, as it is typically the case for a strong pulse (for 
3 X IQi'^W/cm^ and 800 nm wavelength, 



"^quiv 



500 ao 

while a is of the order of lao). However, at 10 times 
higher photon frequency and a factor 100 weaker peak 
intensity which is realistic, e.g., for the FLASH source in 



Hamburg, Xq 



o. Even more problematic is the fact 



that the soft-core potential is harmonic about its min- 
imum with a characteristic frequency (Z/a'^)^/^ which 
could become resonant with the VUV laser frequency. 

To avoid these problems we use a different approxima- 
tive potential which has the correct asymptotic Coulomb 
behavior at large distances but lacks an eigenfrequency 
since it has a non-zero slope at r = 0, 



V{r) = (1 

r 



(6) 



Here, a is chosen in analogy to the potential ^ discussed 
above. Note that U{r^O) = V{r^O) = -Z/a. For rare- 
gas atoms a is of the order of one: a — 1.74 for Xe and 
a = 1.4 for Ar. 

For the photon frequency used here {hcu = 20 eV) there 
is no qualitative difference using a U-shape (soft-core) or 
a V-shape potential. The subsequent considerations also 



do not depend on the approximate or exact form of the of 
the Coulomb potential. Therefore we will use the generic 
form v(i,j), which could be either of the three options: 

{|ri — Tjl^^ exact Coulomb 

[{r,-rjf + a^)y^ U-shape 
(1 - e-l'-'-'-^l/'^) |r, - rj|-i V-shape 

(7) 

Our numerical examples presented here have been ob- 
tained with the V-shape potential. 



Identification of localized and delocalized 
electrons 



At photon energies comparable or less than the ioniza- 
tion potential of a cluster ion {hui < 30 eV), most of the 
photo-electrons remain in the cluster, i.e., quasi-free elec- 
trons are produced. They thermalize quickly, i.e., form 
a plasma. 

We have to determine which among these quasi-free 
electrons travel all over the cluster and visit many 
ions (delocalized quasi-free electrons) and which revolve 
about a single ion, that is, are effectively in excited states 
about an ion (localized quasi- free electrons). To do so, we 
record the revolution angle (/)(t) of each classical electron 
about its closest ion as a function of time. If the electron 
j moves for two revolutions {(f>j^ — 47r) about the same 
ion a we consider it as localized, and the period of its 
motion T^" is then given by 2tt = (pjaiT^"). 

The average period T''" of all localized classical elec- 
trons sets the time scale for a coarse grained dynam- 
ics. This time scale changes slowly in real time t due to 
changes in number and energy of the localized electrons, 
see Fig. [5J With an initial guess for the first averaged pe- 
riod To = 1 fs, we update T after a time t = Ti according 
to the general sequel 



N 

EE 



2n 



.{Ti. 



(8) 



where the actual time interval used is increased by the 
standard deviation cr^+i of the mean, Ti-\-\ — Ti^\^(7i^\, 

The coarse graining of time through the time intervals 
Ti, whose length is defined by the localized electrons, 
plays a crucial role for the description of the entire clus- 
ter dynamics. It provides the natural time scale to in- 
terpolate between the explicit time-dependent dynamics 
of the classical electrons and the time-averaged rate de- 
scription of the bound "quantum" electrons. Over an 
interval < t < to + Ti in time all processes involving 
quantum rates will be considered within a fixed cluster 
environment with properties averaged over the previous 
interval [to — Ti_i,to)- 
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FIG. 2; Time intervals Ti of the averaged localized motion of 
the screening electrons as a function of real time for an Ari47 
cluster with laser parameter as specified in Sect. IIVI 



III. COARSE GRAINED PHOTO-IONIZATION 
INTO THE PLASMA 

We are now prepared to calculate atomic properties in 
the environment of other cluster ions and electrons with 
the understanding that all these processes are for a spe- 
cific time interval Ti as introduced in the previous section. 
The photo-ionization dipole matrix elements for many- 
electron atoms provided within the Hartree-Fock approx- 
imation allow one to determine the photo cross sec- 
tion for ionization of individual occupied orbitals to the 
continuum, see Sect. Ill Al 

To apply these cross sections, we have to approxi- 
mately map the present situation of a cluster ion sur- 
rounded by the localized electrons and other charged 
particles (ions and delocalized electrons) into an effec- 
tive single ion scenario. This requires first to determine 
the electronic energy of an ion with its localized electrons 
and then to construct an energy-equivalent configuration. 

A. Electronic energy of the ions including localized 
electrons 

Averaged over Ti we calculate the number ria of elec- 
trons localized about ion a and their mean energy 

/ 2 \ no 

j=l \ / ]>k=l 

with v{j,k), the interaction potential Eq. ([7]) between 
two particles of unit charge at positions Vj and r^, g^t, 
the charge of the ion a and E{qa) the energy of its bound 
electrons. 

The localized electrons are in excited states of the ion 
a, as shown in Fig. [S] Starting from the energy E* of 
this actual configuration, we include them in the photo- 
ionization process by constructing the equivalent config- 
uration of ion a. In this configuration, we relax all local- 



FIG. 3: Sketch for the construction of the equivalent elec- 
tronic configuration (right panel) from the actual configura- 
tion (left panel) of a cluster ion with three localized electrons 
(red) around a 3s^3p^ (upper row) and 3s^3p^ (lower row) 
configuration. Holes are shown as open circles. 

ized electrons but one onto the last occupied orbital of 
the actual configuration. We put the remaining electron 
on a Rydberg orbit, whose energy is given by the con- 
dition that the actual and the equivalent configuration 
have the same energy E* . 

B. Ionization potential of the equivalent 
configuration 

To find the binding energy of electrons in occupied or- 
bitals in the presence of a Rydberg electron, we assume 
that neither the quantum number n nor the angular mo- 
mentum i of the Rydberg electron change upon release 
of an electron from a deeper orbital. Then the ionization 
energy AE is given by 

AEiC^,n)^E{C^+i,n)-E{C„n), (10) 

where E{Cq, n) is the energy of the valence shell config- 
uration C with charge q and an additional electron in a 
Rydberg orbital n. We have omitted i as an index in 
PH)) since the following discussion does not depend on it. 
As an illustrative example we will use £ = 1. The quan- 
tum number n can take values from no up to oo, with 
no corresponding to the situation where the "Rydberg" 
electron is in the lowest possible state; for Argon and 
£ = 1 it is TiQ = 3. Figure |4] provides energies for six con- 
figurations and n = 3 . . . 9. In five of the configurations 
the valence shell does not contain electron holes. Thus 
they are solely defined by the charge: Cq = q. The case 
6*2 =3s3p*np, also shown in Fig.|4l is an example for an 
exception with a hole in 3s. Obviously we have 

E{Cq, CX3) = E{Cq + i,no) = E{Cq,no) + E,^{Cq) , (11) 

i. e. the asymptotic energy n — > oo of a Rydberg series 
coincides on one hand with the origin of the Rydberg 
series of the next higher charge state, as can be seen in 
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FIG. 4: Total energies (circles) for ions with one of the elec- 
trons in a Rydberg state with quantum number n and angular 
momentum The valence shell configuration Cq is speci- 

fied for each set of energies. Fitted curves (thick lines) accord- 
ing to Eq. . As an example the energy E* for an equiv- 
alent configuration (dotted) as defined in Sect. lIIlXl and the 
corresponding energy Ai?* for ionization (dashed) are shown. 
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FIG. 5: Sketch of the energy balance for the photo-ionization 
of a cluster ion a. The cluster potential is represented by 
a single neighbouring ion. The thick full line indicates the 
full cluster potential, the thin long dashed line represents 
the cluster environment, i.e. the cluster potential without 
the contribution from the ion itself. The dotted line is the 
Coulomb potential for + 1, i.e. the field of the ion if the 
electron would be ionized. The interaction of the bound elec- 
tron with the nucleus is represented by A_B*. The saddle 
point of V(r) — {qa -I- l)u(a:, r) defines the position rbar of the 
barrier. 



Fig.m On the other hand it is equal to the sum of the 
ground-state energy and the ionization potential Eyp for 
an ion with configuration Cq. For finite values of n we 
aT^roximate the energies by a quantum defect formula 

i?(C„n) = i?(C„^) - \ i^^^ , (12) 
2 \n- [iq) 

where in contrast to the usual ionic charge Z we use an 
effective one, 

Zeff = (no - Mg) (2i;(Cg, c^) - 1E{Cq, no))'/' , (13) 

chosen such that the first level (no) of the series agrees 
with the exact value while the quantum defect \iq is fit- 
ted. Hence, Eq. \Vl\ is very accurate at intermediate n, 
where we need it. When fitting the curves Eq. p2|) to 
the calculated energies, cf. Fig.[51 we found the quantum 
defects [iq to be almost independent of q. This allows us 
in the calculation of energy differences Eq. (llOp to elim- 
inate the term containing n and in Eq. (|12p. We get 
for a configuration Cq with an initial energy E* 

AE*{Cq) = [E{Cq+i, oo) ^E*]-b [E{Cq, oo) - E*] , 

(14) 

with 

6 ._ EipjCq+l) ^ ^^^^ 
Eip{Cq) 

the ratio of the ionization potentials, cf. Eq. (fTT|) . Thus 
we have obtained an expression for the energy necessary 
to ionize an electron from a valence shell in presence of 



a Rydberg electron. This expression does not depend on 
the actual quantum numbers n and £, but only on the 
energy E*. 

C. The condition for over-barrier inner-ionization 

Although we know now with Eq. (jl4p the ionization 
potential for a screened isolated ion, we have to position 
the ion in the cluster environment in order to decide, if 
photo-absorption leads to photo-excitation within the ion 
a or to inner-ionization above the lowest barrier on the 
way to a neighbouring ion. The energy balance for the 
photo electron which decides between these two options 
is 

Ei + Vi+Lu = Ef + Vf (16) 

with local contributions from the ion to which the elec- 
tron is bound, (see Fig. [5]) 

E,^~AE*, =^ki„- (9a + l)f(a,bar), (17) 

and contributions from the background of charges in the 
cluster, 

Vi^ - ^ qjv{j, a), V/ = - ^ qjv{j, bar) , (18) 

where the index j runs over the delocalized electrons and 
all ions but ion a. As introduced earlier, v(l,m) is the 
interaction between two Coulomb particles at positions 
r/ and r^. Hence, Vi refers to the potential energy of the 
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FIG. 6; The number of localized quasi-free electrons nioc and 
the corresponding energy E* according to Eq. © for the cen- 
tral ion and a surface ion of an Ari47 cluster from the micro- 
scopic calculation in Sect. II Vl 



electron under consideration and located at the position 
of its mother ion a due to interaction with particles of 
charge qj at rj. Likewise, Vf is the potential energy 
of the same electron at the potential barrier rbar due to 
the interaction with the same charged particles as before. 
The energy balance Eq. is taken with respect to the 
location of the ion Tq, and the location rbar of the lowest 
potential barrier near the ion. AE"^ in Fig. El is defined 
by putting ii^kin = in Eq. ^7^, i.e., 

A£:'=*^ = AE* - (fc + l)v{a, bar) + Vf-V^. (19) 

Figures [6] and [7] give an overview of the coarse grained 
variables during the laser pulse for an atom at the cen- 
ter and one at the surface of the Ari47 cluster. Figure [5] 
shows first of all an overview of their evolution during the 
whole interaction with the laser. Starting at the ground 
state of the neutral Ar, each absorbed photon leads to a 
rising step in the total energy. Note that for an ion in 
the cluster the electron has to be excited only above the 
lowest barrier. Moreover, the energy E* of the equiv- 
alent configuration takes merely the localized electrons 
into consideration and not the newly ionized one. There- 
fore, each ionization event leads to jumps of E* higher 
than the energy of the photon, cf. Fig. [6l 

The electron localization is equivalent to a relaxation of 
the system and lowers therefore E*. The fiat regions ob- 
servable for both ions correspond to the case where there 



-14.20 



o central atom 
X surface atom 




143 144 145 146 



147 
t [fs] 



148 149 



150 



151 



FIG. 7: Zoom into Fig. [6] for 143 < t < 151 fs. The upper 
panel shows the number of localized electrons (o for the cen- 
tral atom, where g = 6 and x for the surface atom, q — 4). 
The lower panel shows the total energy E* of the ions and 
their localized electrons. The full lines show here the average 
of E* over the time intervals Ti, as introduced in Section FlI CI 



are no localized electrons, when the total energy of the 
ion is given solely by the bound, "quantum" electrons. 
The smaller final charge of the surface atom is a con- 
sequence of the cluster expansion. The surface expands 
much faster than the core, leading to higher interionic 
barriers and an early suppression of the inner-ionization. 
A detailed view of the evolution of the two atoms is shown 
in Fig.[7]for the time interval from t = 143 to 151 fs. The 
coarse graining is symbolized by the full lines, showing 
the total energy E'* of the two atoms averaged over the 
time intervals Ti, as described in Sect. Ill CI 



D. The effective cross section for inner 
photo- ionization 

Finally, we are in a position to adopt the photo- 
ionization cross section as formulated in Sect. Ill Al for 
an isolated ion to the situation of an ion in the cluster. 
Here, we take the lowest potential barrier to a neighbor- 
ing ion as effective ionization threshold. Therefore, the 
actual cross section as shown in Fig. [T] can vary between 
minimum and maximum possible values of potential bar- 
riers. The interval is indicated for small photon energies 
in Fig. [1] with the additional dotted line. 

The electrons available for photo-ionization are only 
the tightly bound ones from the actual configuration, 
while the matrix elements entering the expression for the 
cross section Eq. ^ take into account the screening of 
the electrons as provided by the equivalent configuration 
cr ~ cr ((7oqv)- Hence, the multiplicity has to be taken from 
the actual configuration w^ct to arrive at the screened 
photo cross section 



'■(g) ^ Zf^cr(q^q^) 



(20) 



cqv 
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IV. DYNAMICS OF AR147 UNDER AN 
INTENSE VUV LASER PULSE 

We will illustrate the theoretical framework introduced 
above with the dynamics of Ari47 exposed to a 100 fs 
VUV laser pulse with huj = 20 eV (i. e. wavelength 
A = 62 nm) and an intensity of 7 x 10"^^ W/cm^. Figure[8] 
shows the main quantities characterizing the response of 
the cluster to the VUV pulse, namely the energy absorp- 
tion, the ionization degree, the temperature of the plasma 
and the cluster explosion. The latter is characterized by 
the interionic distance i?avgi defined as 



1 

A — > minflrfl 

\ a— 1 



(21) 



In the first part of the pulse, until approximatively t — 
80 fs, the absorption is dominated by photo- ionization. 
After that, inverse bremsstrahlung (IBS) sets in and 
photo-ionization starts to saturate. We can easily disen- 




400 



FIG. 8: Explosion of Ari47 exposed to a 100 fs long VUV 
laser pulse with intensity 7 = 7 x 10^^ W/cm^ and photon 
energy hu; = 20 eV. (a) full-line: total absorbed energy per 
atom; dashed-line: energy absorbed due to photo-ionization; 
(b) full-line: total number of ionized electrons per atom; o 
quasi-free electrons; o localized electrons; (c) temperature of 
the quasi-free electrons; (d) expansion: average interionic dis- 
tance, see Eq. (HTJ. 



tangle the two contributions to the total absorbed energy 
since we treat photo-ionization via rates, see Sect. Ill A[ 
and IBS through classical propagation, see Sect. Ill Bl 

Note that until this time, as can be seen in panel 
(b) of Fig. [H almost all electrons (ritotai) are trapped 
('T-quasi-ficc) iusidc the cluster due to their low kinetic 
energy. They start to leave the cluster with the onset of 
the IBS heating. This disagrees with the common inter- 
pretation that, when the photon frequency is larger than 
the first atomic ionization potential, all cluster atoms 
are singly ionized by one-photon absorption. Rather, the 
space charge built up even in a relatively small cluster 
such as Ari47 allows the electrons to become quasi-hee 
only. Hence, the "ionization into the cluster" , as de- 
scribed in Sect. 1111) starts very early in the pulse. These 
electrons thermalize very quickly and obey a Maxwell- 
Boltzmann velocity distribution with a temperature T 
which is shown as function of time in Fig. [8t. 

Figure [8)d also shows the average number of localized 
electrons, as they have been introduced in Sect. Ill CI 
They reach a maximum (t k, 90 fs) when the photo- 
ionization becomes unlikely. The increase of the temper- 
ature of the electron plasma (see Fig.[5J:) from this point 
on favors a decrease of the number of localized electrons. 
At the point {t k. 120 fs) where the cluster expansion (see 
Fig. IHli) startes, the temperature of the electrons plasma 
decreases, despite of the continuing energy absorption 
and the localization increases again to an average of two 
electrons per atom. 
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FIG. 9: The ion charge distribution (left panel) and the ki- 
netic energy distribution of the released electrons (right) af- 
ter the cluster explosion for the same pulse as in Fig. jS] The 
propagation time was t = 6.4 ps. The electron signal can 
be exponentially fitted by exp(— iJkin/So) with Eq = 5.4 eV 
(dashed line). 

This transient complex behavior of electrons and ions 
in the cluster may be accessible experimentally soon us- 
ing attosecond pulse probe techniques, for a proposal see 
[23'|. For the time being, we present in Fig. [9] more con- 
ventional observables measurable in the experiment, such 
as the final charge distribution of ions and the kinetic- 
energy distribution of electrons for the dynamics of Ari47. 
They have been obtained by propagation up to i = 6.4 ps, 
i. e., much longer than shown in Fig. [S] The dominating 
fragment is Ar^+ despite the fact that about five electrons 
(Fig.[5}3) per atom have been photo-ionized. The kinetic- 
energy distribution of the electrons can be fitted with 
an exponential decay exp(— i^km/^-o) with = 5.4 eV, 
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thus emphasizing the thermal origin of the electrons. The 
temperatm'e of the plasma has similar values, see Fig.[51:. 

V. SUMMARY 

We have described a quantum-classical hybrid ap- 
proach to follow the time-dependent dynamics of rare- 
gas clusters exposed to intense short laser pulses. Spe- 
cial attention has been paid to incorporate the screen- 
ing of cluster ions by the electron plasma, formed by 
quasi-free electrons which have been ionized from their 
mother ions but cannot leave the cluster as a whole due 



to the strong background charge. The mean time scale 
of these localized quasi-free electrons is determined and 
provides the link between the microscopic dynamics of all 
charged particles and the quantum dynamics of photoion- 
ization which is described by ionization rates adopted to 
the screening in the cluster environment. 

Hence, this approach is especially well suited to tackle 
interaction of clusters with light from VUV and X-FEL 
laser sources. As an illustrative example we have dis- 
cussed the dynamics of Ari47 exposed to a 100 fs laser 
pulse of / = 7 X 10^^ W/cm^ intensity and ftw = 20 eV 
photon energy. 
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